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Abstract 

Motivated by ATLAS diboson excess around 2 TeV, we investigate a phe¬ 
nomenology of spin-1 resonances in a model where electroweak sector in the SM is 
weakly coupled to strong dynamics. The spin-1 resonances, W and Z', are intro¬ 
duced as effective degrees of freedom of the dynamical sector. We explore several 
theoretical constraints by investigating the scalar potential of the model as well as 
the current bounds from the LHC and precision measurements. It is found that the 
main decay modes are V —>■ VV and W —)• V/i, and the V width is narrow enough 
so that the ATLAS diboson excess can be explained. In order to investigate future 
prospects, we also perform collider simulations at t/s = 13 TeV LHC, and obtain 
a model independent expected exclusion limit for a{pp —>■ W —)• WZ —)• JJ). We 
hnd a parameter space where the diboson excess can be explained, and are within 
a reach of the LHC at f dtC = 10 fb“^ and ^/s = 13 TeV. 
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1 Introduction 

The ATLAS experiment recently reported an excess of the events in the search for the 
diboson resonance, in the pp —?• (hPlT, hPZ, and ZZ') —)■ JJ channels, where J is a fat- 
jet formed by boosted IT or Z boson [1]. The largest local signihcance is 3.4 a around 

2 TeV in the WZ channel, and the global signihcance is 2.5 a. The CMS experiment also 
studied the same channels. The largest deviation they found is 1.4 cr at ~ 1.9 TeV [2]. 
Although we cannot conclude that there is a new particle with the mass around 2 TeV 
from this data, it is worthwhile to consider models which can explain this excess, and 
many papers have already appeared discussing interpretations of the excess [3-18]. As 
discussed in these references, a simple candidate is a spin-1 particle. 
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New vector resonances often appear in the models with dynamical symmetry break¬ 
ing. Such spin-1 resonances appear in the composite Higgs scenario [19-22] with the 
dynamics at TeV scale to account for the naturalness problem. Since the models are 
based on the non-Linear sigma models, the effective theory involves many operators 
whose coefficients are unknown. They possibly affect to the couplings of the new spin-1 
particles to the standard model (SM) particles, and thus there is uncertainty in the 
prediction of the W couplings. 

Another way to include spin-1 particles is to extend the electroweak gauge symme¬ 
try. We can easily introduce new spin-1 particles in renormalizable manner. In that 
case, the models are calculable and we can avoid the operators whose coefficients are 
unknown, in contrast to the models based on the non-linear sigma models. Besides, 
some renormalizable models with extended gauge sector can be regarded as the low 
energy effective theory of UV theory with some dynamics. 

Such renormalizable models have been discussed in the context of the left-right (LR) 
symmetric model [8,17,18] and the leptophobic G221 model [11]. These models contain 
the right-handed SM fermions which are not singlet under the new gauge symmetry. In 
such case, the couplings of the SM fermions to the new gauge boson are not suppressed, 
and the new gauge bosons mainly decay into the SM fermions. 

It is also possible to use linear sigma model, instead of non-linear sigma models, for 
models emerged from the dynamics at TeV scale. An example was proposed in Ref. [23]. 
This model, called the partially composite standard model, has three Higgs helds. Two 
of them are regarded as effective degrees of freedom below the dynamical scale. The 
other one is an elementary held. Spin-1 resonances are introduced as new gauge bosons 
a la Hidden local symmetry [24-27] . A feature of the model is that the SM fermion are 
singlet under the new gauge symmetry, and thus all the fermion couplings to the new 
gauge bosons are suppressed by the mixing angle in the gauge sector. As a result, the 
new gauge bosons mainly decay into the SM gauge bosons. This is an important feature 
of this model. 

In this paper, we investigate the possibility to explain the diboson excess by the 
partially composite standard model, and also the future prospects of W and Z' bosons 
searches at the LHC Run-2, where ^J~s = 13 TeV. We perform a comprehensive study 
to hnd the parameter space which has not been excluded from current experimental 
data. The constraints on the model parameters come from the LHC searches and the 
electroweak precision measurements. We also require theoretical constraints such as 
perturbativity condition, bounded below condition, global minimum vacuum condition, 
and stability condition of the scalar potential. We hnd a parameter space where the 
diboson excess can be explained. 

We further investigate a model-independent sensitivity at the LHC Run-2 by gen¬ 
erating both signal pp —)■ W WZ and dijet background events, and performing 
detector simulations. 

We organize the rest of this paper as follows. We briehy review the partially compos¬ 
ite standard model in Sec. 2. In Sec. 3, we show the constraints to the model, and hnd 
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Figure 1. The moose diagram of this setup: the circles represent the gauge groups, and the 
thick lines that connect two circles are the Higgs fields. The Higgs fields Hi, H 2 , and SU(2)i 
gauge group can be regarded as the operators originated from the dynamical sector. 


that there are parameter regions where are consistent with the ATLAS excess. In Sec. 4, 
we perform the collider simulations for the signal and the background, and obtain the 
sensitivity at the LHC Run-2. In Sec. 5, we investigate the future prospects of the spin-1 
resonances search using our simulation results. Section 6 is devoted for conclusion. 

2 The partially composite standard model 

2.1 The model setup 

In the partially composite standard model, the gauge symmetry of the electroweak sector 
is SU( 2 )oxSU( 2 )ixU( 1 ) 2 , and three Higgs helds (Hi, H 2 , H 3 ) are introduced for the 
symmetry breaking, SU(2)o xSU(2)ixU(l)2 —)• U(1)^m- We denote the gauge couplings 
go, gi, and g 2 , respectively. The three gauge couplings are related to the QED coupling 
as 
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( 1 ) 


Here, we assume that the SU(2)i gauge symmetry belongs to a dynamical sector, and 
91 3> go, 92 - Under this assumption, go and g 2 are approximately gw and gy which 
are the gauge couplings of SU(2)i and U(l)y, respectively. We regard the gauge held 
associated with SU(2)i as the vector resonance originated from unknown dynamics above 
TeV scale. This implies that helds transformed under the SU(2)i gauge symmetry also 
belong to the dynamical sector. We take Hi and H 2 as such helds, and regard H 3 as 
an elementary held. All the fermions are also elementary, and they are singlet under 
SU(2)i. We schematically show the model structure in the moose notation [28] in Fig. 1, 
and also summarize the held contents and their charge assignments in Table 1. 
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Table 1. The charge assignment of the partially composite standard model. Only Hi and H 2 
are the representations of the SU(2)i gauge symmetry. 


Fields 

SU(2)o 

U(l )2 

SU(3), 

SU(2)i 

Hi 

2 

0 

1 

2 

H 2 

1 

1/2 

1 

2 

Hs 

2 

1/2 

1 

1 

Q 

2 

1/6 

3 

1 

Ur 

1 

2/3 

3 

1 

Ar 

1 

-1/3 

3 

1 

L 

2 

-1/2 

1 

1 


1 

-1 

1 

1 


The scalar potential is given as"^^ 

V{Hi, H2, H^) = filtriHiHl) + 

+ KtT{HiH2Hl) 

+Xi{tT{HiHl)f + X2{tTiH2Hl)f + XsitiiHsHl))^ 

+ Xi2tT{HiHl)tT{H2Hl) + X23tT{H2Hl)tT{HsHl) 
+X3MH3hImh,hI). ( 2 ) 

Here all the Higgs helds are represented by two-by-two matrices, and they are real, 
namely 


eH*e = -Hi, where e = . (3) 

All parameters in the Higgs potential are also real. We assnme that all the vacnnm 
expectation valnes (VEVs) of the Higgs helds are diagonal, real and positive to realize 
desired electroweak symmetry breaking. The Higgs helds are expanded aronnd their 
VEVs, vi, V 2 and V 3 , 


«. = I + ^ (fc. + , (4) 

where r“ is the Panli matrices, and hi, 7 r“ are the fonr real scalar component helds. The 
covariant derivatives of the Higgs helds are given as 

D,Hi = d.Hi+tgu^WS.Hi-tgiHi^W^^, (5) 

We omit the term in'tv{HiH 2 H\H) in this paper because this term can be eliminated by the field 
redefinition [23]. 
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( 6 ) 


D^H2 = d^H^ + igJ-W-^H^-ig^Hj-B^, 
D^Hs = d^H. + tgo^WS^H.-tg^H.^B^, 


(7) 


By calculating the muon life time in this model at the tree level, we hnd the relation 
between the Fermi constant and the VEVs in this model as 


+ T-r = 


— -I- — 
..2 \ .,2 


-1 


( 8 ) 


where v ~ 246 GeV. For the later convenience, we introduce a new parameter r. 


r = — 

~ Vl' 


Thus Vl and V 2 are expressed by r, U 3 , and u, 

«? = (! + T-^)(V^ - vl), vl = (1 + - vl). 


(9) 


( 10 ) 


There are twelve scalars in this model, and six of them are eaten by the gauge bosons. 
Thus this model has six physical scalars: three CP-even Higgs bosons (h, H, H'), one 
CP-odd Higgs boson (H), and two charged Higgs bosons {H^). We identify h as the 
SM-like 125 GeV Higgs bosons. The masses of the CP-odd and the charged Higgs bosons 
are the same at the tree level and given by 


2 2 I K 1 + 2 

rriA = mTT± = —- v . 

^ " 4 U 3 r 


(11) 


The mass eigenstates of the CP-even Higgs bosons are related to the gauge eigenstates 
through the mixing angles 6 ^ 1 , 6 ^ 2 , and 6 ^ as follows. 


( H' \ / S1S2 - C1C2S3 -S1C2 - C1S2S3 C1C3 

if = -C1S2 - S1C2S3 C1C2 - S1S2S3 S1C3 

\ h J \ C2C3 S2C3 S3 

where Sj (c*) stands for sin 6 *j (cos 6 *j) for i = I, 2, 3. 

The Yukawa interactions are given as 


hi 

h 2 

hs 


C 


Yukawa 


Q'H, 






'R 


+ H.C., 


( 12 ) 


(13) 


where i and j are the generation indices. We introduce a parameter Kp which is the 
ratio of the couplings between the lightest CP-even Higgs boson and the fermions to its 
SM value. 


Up 


- 9hff 
rrif/v 


V 



(14) 
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Since jssl < 1 and v/v^ > 1, Kp can be larger than one. We will discuss the viable range 
of Ki? in the next subsection. The fermion masses are given as 

^3 ^3 V 

( 15 ) 

and the Yukawa couplings are enhanced by v/v^, compared to their SM values. Large 
Yukawa couplings could make the Higgs potential unstable above the electroweak scale. 
We discuss this point in Sec. 3.2. 

In addition to the SM gauge bosons, we have extra three vector bosons, W'^ and 
Z'. In the Qi 3> go regime, the mass eigenvalues of the gauge bosons are given as 


m| ~ 

m|/ ~ 


1 

4 

1 

4 

1 

4 

1 

4 


2 2 
9o^ 



9l 1 \ 




9l{vl + vl) 

{ 9 l + 9 iy 

9l{vl + vl) 


fl + gl ^ ^ 

A {9l-9y? 1 \ 

V 9l{9l + 9D ii + yV' 

A , 9l + 9y 1 \ 

V 9i {i + yy' 


(16) 

(17) 

(18) 
(19) 


We hnd mw' — except in the large r regime. We need to hnd the relation between 
the gauge eigenstates and the mass eigenstates to evaluate the couplings, and they are 
given as 





2(1 + r2)2 gf 

1 9o 


1 + 


r^gi 


K + 



e n e n e n 


2(1 + r2)2 gl 



9o 

Zfj, — Cw 


92 

-2 ^2 


■ Sw 


_ 1 - 2rH^w 

2 (l + r 2 ) 2^2 

' i/o 


+ Cw 



2(1 + r 2)2 g^ 




Z'a-- 

M 1 I 


1 9o 




f _ l + rHl. 9o \ ...3 _ 9o 

W on I »,2^2 „2 n I »,2^ „ 


T;r/3 


( 20 ) 

( 21 ) 

( 22 ) 

( 23 ) 

( 24 ) 


where cw = mwlmz, sw = \/l — c^, and tw = sw/cw- The typical size of the mixing 
angles is O{go/gi), and the gauge hied ITo (lYi) is the main component of the mass 
eigenstate W iW'). It is worth noting that the mixing angles for W and Z' become the 
same in the small r regime, which means the custodial symmetry is enhanced. 
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The approximate expressions for some couplings of the W and Z' to the SM particles 
are given as 


9w'ff 

mw 

h- 

nil 

g^ff 

mw' ' 

r 

yA j. 

gw'wz 

mw 

/.- 

ni 1 r 

9\^z 

mw' ' 

n2 1 + r2 

gwwz' 

mw 


ni 1 r 

9v^z 

mw' ' 

n2 cvi/1 + 


(25) 

(26) 
(27) 


Compare to the benchmark model (sequential standard model (SSM) [29-31]) used by 
the ATLAS/CMS, the W couplings to the SM fermions are smaller. All couplings have 
a suppression factor of {mw/^w')\/^ — because the W boson couples to the 

elementary fermion through Wq-W\ mixing, so that the width is narrow. Since gw'ff is 
proportional to the production cross section of the Drell-Yan process is proportional 
to r“^. The W boson could not be produced in the large r region. In Sec. 3, we will hnd 
a parameter space in small r regions where a signal rate is consistent with the ATLAS 
diboson excess. Some numerical results are given in Sec. 3.3 

The gauge boson couplings to the scalars are also important in our analysis. We give 
their approximated expressions here. Due to the mixing between the two SU(2) gauge 
eigenstates, the WWh and ZZh couplings differ from the SM values. We denote these 
coupling ratios to the SM values by 


_ QhWW _ dhZZ 

2 m^/n’ 2 m|/n’ 

and their approximated formulae in the gi 3> go,g 2 regime are 


(28) 


KW - I^Z - 


(l + r2)3/2 



C 2 C 3 + 


1 

(l + r2)3/2 



vl v-i 

4 s 2 C 3 H-S3. 

V 


The couplings relevant to the W jZ' decay are 


(29) 


2mwmw' 

gww'h ^ gzz'h ^- 

V 


2mwmw' 

gww'H - gzz'H — - 

V 


^2)3/2 ^2^3 + ('X + ^2)3/2 '^2^3 ^ 


1 V2, 



5 


(1 + j-2^3/2 


(-C1S2 - S1C2S3) 


r 

+ (l + r2)3/2 


(CiC2 - S1S2S3) 


1 V 3 
rrri^, v 
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2mw'mw' 

gww'H' — QZZ'H' — - 

V 


(l + r2)3/2 


(SiS2 ~ C1C2S3) 


(1 + r2)3/2 


[ — S1C2 — C1S2S3) 


1 V 3 

rrn^, v 




(32) 


2.2 Model parameters 

In the electroweak sector of the model, there are 13 real parameters, 

hi) hL hs) ^) -^1) -^2) A3, Ai2, A23, A31, 5^0) 91 , 9 i - (33) 

It is convenient to use a different set of the parameters instead of these parameters. We 
use the following 13 parameters to £x the parameters in the electroweak sector, 

r. Vs, V, a, mz, ruz', ruh, tuh', mu, mA, Kp, kz, Qww'h'- (34) 


Here, we use r, ^ 3 , v instead of the three /i parameters, k is fixed by the charged Higgs 
mass or the CP-odd Higgs mass. The six A’s have the same information as the three 
CP-even Higgs masses and their mixing angles (m^,, mp, thh', Oi, 62 , O 3 ). We can use 
mz, mz', and a(= e^/dvr) instead of the gauge couplings. In addition, we can replace 
( 6 * 1 , 6 * 2 , 6 * 3 ) with (up, Kz, gww'H’)- 

The values of the four parameters, v, mz, a, and are already known very pre¬ 
cisely. We take = 125 GeV. We find np is severely constrained close to 1 in the 
mA S> mh regime. When the heavy Higgs masses are universal, mA = mp' = mp, A 3 is 
simply expressed as 


A3(/i = mz') 


K^p m\ 1 — K^p m\ 
^ 2 


(35) 


Thus, A 3 is very large except for Kp =1. Similarly, the coupling ratio kz is also severely 
constrained close to its maximal value in mA ^ mh regime. Typically, the allowed value 
of is 1 — 0{1)%. The detailed description is given in Appendix A. In the following 
discussion, we take Kp ^ 1, mA = mp = mp/ = 0(1) TeV, gwwH' = 0, and we always 
choose Kz to be its maximal value. For most of our numerical analysis, kz is set to 
0.95-1.00. 


3 Phenomenology of spin-1 resonances 

In this section, we discuss the properties of the W and Z' bosons such as the production 
cross sections and the decay branching ratios. And we also discuss both the theoretical 
and experimental constraints. In the following, we show some formulae with approxima¬ 
tion, which help to understand the parameter dependence. However, we use the exact 
formulae in our numerical calculations. 









Figure 2. The production cross sections of the extra vector bosons, W and Z'. We take 
ns = 200 GeV, r = 0.13. 


3.1 Properties of the extra vector bosons 

The main production mode of W and Z' in pp collisions is the Drell-Yan process. They 
are proportional to as we can see from Eq. (25)."^^ We consider small r in order 
to make the cross section large. In Fig. 2, we show the production cross sections for 
W (IT'+ + W'~) and Z'. Here we take = 200 GeV, r = 0.13 and Kp = 1.00. We use 
the CTEQ6L parton distribution functions [32]. The production cross section of W 
is approximately twice that of Z' for small r because the custodial SU(2) symmetry is 
recovered in the region. 

The partial decay widths of W and Z' into the SM particles are given as 


V{W' WZ) 
T{W' Wh) 

r(w' ^ //) 

T{Z' WW) 
T{Z' Zh) 

r(z' ^ //) 


1 


y2 


dSvr 

v2 

(1 + 

y 

1 


y2 


dSvr 

v2 

(1 + 

r 

Nc 


e2 1 1 

1 

00 

1 

mw' 



1 

TTV^, 

y2 


dSvr 

v2 

(1 + 


1 

Tfl^, 

y2 


dSvr 

v2 

(1 + r2 


Nc 

rri^ 

1 / 

f 

1 

247r 

Tn\Y' 

s^r2 \ 



1-^ 


{-rC2C3 + S2C3f , 




i-d 


{-rC2C3 + S2C3f , 


^^For more details on the formulae, see Appendix B. 


(36) 

(37) 

(38) 

(39) 

(40) 
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X 


(41) 





+ 



1 


where T| = 1/2 (—1/2) for the up-type (down-type) fermions, Qf is the electric charge 
of the fermion, and Nc is the color factor, Nc = 3. Note that the bosonic channels 
are dominant among their decay channels due to an enhancement factor from the wave 
function of the longitudinally polarized gauge bosons in the final states. Therefore the 
dominant decay modes of W {Z') are W ^ WZ and W ^ Wh {Z' ^ WW and 
Z' Zh). 

There are also decay modes with a heavy scalar in the final states. For example, 
W WH and W —)■ hH exist and their approximated formulae are 


T{W' WH) ( 


Qww'h 


dSvr \2m-\vmwi/v J v‘^ \ m 


m 


W' 


m 


w 


+ 1 


‘'W' 



m 


H 


m 


w 


(42) 


T{W' hH) ~ 


m 


w 


487r (1 -|- r^)^ 


(-rC2C3 -1- S2Czf ( 1 


m 


H+ 


m 


W' 


(43) 


where gwwH is given in Eq. (31), and gww h// v) is 0{1) in large regions of 
the parameter space. These decay widths are comparable to those of the SM final states. 
The decay channels including only the heavy states, such as W —)■ also have the 

same feature. Once these modes are open, they would be dominant decay modes. 

In Fig. 3, we show the total widths and the branching ratios of W and Z'. Here we 
take U 3 = 200 GeV, r = 0.13 and = 1.00, = rnn' = mn = 2 TeV. We also take 

gww'H' = 0 as we mentioned in Sec. 2.2, thus W WH' and Z' —)■ ZH' are absent 
in the hgure. We hnd that the dominant decay channels are V —)■ VV and V ^ Vh 
for large my'- However, the decay branching ratio reduces once V —>■ HX decay modes 
are open. 

In Fig. 3, we find Br(iy' ^ WZ) = W{W' Wh) ~ 40 % at mw = 2 TeV. The 
relation is easily understood by the equivalence theorem, Br(hF' —?• WpZp) =Br(kF' —)■ 
tism'^sm) in the heavy W' mass limit, where t^sm is the SM Nambu-Goldstone boson. 
In the SM limit, Br(kF' —>■ ttsm'^sm) is equivalent to Br(IF' —)■ TTsuh). Thus, Br(IF' —)■ 
WZ) ~Br(IF' Wh) is realized. In addition, when one takes gwwH' = 0, the heavy 
Higgses H, A, H^ form a multiplet, and Br(IF' —)■ ZH^) ~Br(IT' WA) ~Br(IV' —)■ 
WH) is realized. 

We also hnd the total widths are narrow, namely F/m ~ 1-5 %, because the reso¬ 
nances decay into the SM particles with suppressed couplings due to the mixing between 
the elementary and the composite sectors. Another remark is that the decay properties 
of W' and Z' are similar due to the enhanced SU(2) custodial symmetry in small r 
regime. 

the large r region, the branching ratio Z' ff and the production cross section of Z' become 
larger (see Eqs. (41, 72)), while the production cross section of W is suppressed. 


10 














mw[GeV] mz.[GeV] 



Figure 3. The total widths and the branching ratios of the extra vector bosons, W' and Z'. We 
take U 3 = 200 GeV, r = 0.13, kf = 1.00, and niA = rnu' = fnH = 2 TeV. Here, 2jets means 
Br(VF' —)• ud)+'Qi{W' —)• sc) or Br(Z' —)■ uu)+Br(Z' —)• dd)+'Bi{Z' —)> ss)+Br(Z' —>• cc), 

means Br(VF' —)• eVe) ( =Br(VF' —=Br(VF' — tUt)), ii means Br(Z' —>• ee) ( 
=Br(Z' —)• an) =Br(Z' —>■ rr)), va means Br(Z' —)> r'et'e)+Br(Z' —)• z^uJ/u)+Br(Z' —)• 
and means Br(lF' ^ W+R-) =Bi{W' W-R^). 
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Figure 4. The total widths and the branching ratios of W' and Z' for niA = "mH' = innH 
TeV. Other parameter choice and the notations are the same as in Fig. 3. 
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In Fig. 4, we take different scalar masses, niA = mn' = mn = 1 TeV. The other 
parameter choice is the same as in Fig. 3. The decay channels to the non-SM particles, 
namely Br(hF' ^ H+H), Br(hF' ^ H+A), Br(Z' HA), and Br(Z' ^ are 

open. We took the parameter so that the channels to H' is absent, Qww'h' = 0. 

As mentioned above, these channels are comparable to the decay modes to the SM 
particles. As a result, the decay modes V —)■ VV and V ^ Vh searched at the ATLAS 
and the CMS experiments are suppressed compared to Fig. 3. This implies that the 
decay channels with heavy scalars should not be open if we try to explain the excess at 
the ATLAS experiment. Hereafter, we consider the situation that the heavy scalars are 
as heavy as the extra gauge bosons. 

3.2 Constraints on the model 

In this subsection, we show theoretical and experimental constraints on the model pa¬ 
rameters. In order to perform a reliable perturbative calculation, we demand pertur- 
bativity condition, bounded below condition, global minimum vacuum condition, and 
stability condition. In addition, we take into account the LHC bounds and the elec- 
troweak precision measurements. 

3.2.1 Theoretical constraints 

Perturbativity condition for the gauge coupling gi 

We require that all the absolute values of the gauge couplings and the Higgs quartic 
couplings are smaller than dvr and (dvr)^ in order to keep the reliability of our analysis 
based on the perturbative calculation. 

For the gauge couplings, go and g 2 are almost the same value as the SM gauge 
couplings, but gi can be very large. Since SU(2)i is asymptotic free in our setup, gi 
becomes smaller at high energy due to the quantum effects. Thus the maximum value 
of the gi is given at Z' mass scale, and we require \gi{fi = mz')\ < dvr. The regions 
where \gi{mz')\ > dvr are hlled with yellow in Figs. 5 and 6. 

Perturbativity condition for the Higgs quartic couplings 

The scalar quartic couplings are large when the mass differences between the SM-like 
Higgs and the other heavy scalar bosons are large. Since we take mh -C = mn = 
tuh' ~ C>(1) TeV in our analysis, the quartic couplings tend to be large. In addition, 
due to the renormalization group effects, they can be even larger at high energy. The 
renormalization group equations for this model are given in Appendix C. We dehne a 
cutoff scale A by 

|Ai(/r = A)| = (47r)% (44) 

and require |Aj(/i < A)| < (dvr)^. 
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This constraint highly depends on onr choice of A. We reqnire that A shonld be 
signihcantly higher than a few TeV, otherwise we have to take acconnt of interaction 
terms from higher dimensional operators whose coefficients are nnknown, and some nn- 
certainty is introdnced to onr analysis. For example, there are operators which modifies 
the W conpling to the SM fermions snch as {c/K^)Qi'^^{HiiD^H\)Q. This operator 
brings nnknown parameter c, and thns brings nncertainty to onr calcnlations snch as the 
prodnction cross section of W. Typically, snch higher dimensional operators with c ~ 1 
bring 1% (10%) nncertainty if A = 100TeV(10TeV). To avoid snch nncertainty from 
nnknown parameters, we restrict onrselves for the case A > 100 TeV. The parameter re¬ 
gions where |A(A)| > (dvr)^ are filled with the lighter (darker) gray for A = 100(10) TeV 
in Figs. 5 and 6 for reference. 

Bounded below condition 

Here we consider the conditions that the Higgs potential at the tree level is bonnded 
below. For the pnrpose, it is enongh to check that the potential valne at the large field 
valnes, and thns we consider only the qnartic terms in the potential. We rewrite the 
qnartic terms as 

Hquartic = (Ain^ A 2 n 2 AaTig Ai 2 nin 2 X23n2n3 + Asinsni) , (45) 

with 

R^Ui = tT{HiHj) = ^(h^ + «), (46) 

where rij satisfies 0 < n-i 2,3 < 1, and nl + 712 + = 1. In order to avoid rnn-away 

vacna, we demand the following conditions for the Higgs qnartic conplings at m^/ scale 
(Aj(/r = mz')), 


Min [Ain^ A 2 n 2 Asrig Ai 2 nin 2 A23?^2?^3 + A 3 in 3 ni] > 0, for 0 < ^ 1 , 2,3 < 1- (47) 

The parameter regions where this condition is not satisfied are filled with cyan in Figs. 5 

and 6 . Especially, one can solve the above ineqnality analytically in specific directions 
as 


Ai > 0 for 77-2 = 77,3 = 0, A2 > 0 for ni = = 0, 

A3 > 0 for TTi = 712 = 0, A23 -l- 2\/A^A3 > 0 for rii = 0, 

A31 -|- 2 a/A3A1 > 0 for 712 = 0, Ai2 -I- 2a/AiA2 > 0 for 713 = 0. (48) 

Global minimum vacuum condition 

We demand the electroweak vacnnm to be a global minimnm of the Higgs potential at 
/i = mz'- The parameter regions where this condition is not satisfied are filled with 
green in Figs. 5 and 6. 
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Stability condition 

Among the quartic couplings, A 3 can be very small when np is very close to 1, see 
Eq. (35), and especially it takes the same value as the quartic coupling in the SM for 
Kp = 1- In that case, A 3 can be negative at a high energy scale due to the contribution 
from the Yukawa interaction to the renormalization group equations, and the Higgs 
potential becomes unstable. The VEV giving masses to the fermions are U 3 , see Eq. (15). 
Thus the Yukawa coupling in our setup is larger than the coupling in the SM by v/v^, 
and the Higgs potential can become unstable at a few TeV scale for the small U 3 region. 
We dehne the scale A at which A 3 becomes negative, 

A 3 (/i = A) = 0 , (49) 

and we demand A > 100 TeV, as we demand for the perturbativity condition. We £11 
the regions where this condition is not satisfied with magenta in Figs. 5 and 6 . 

This bound is conservative because we do not allow a meta-stable vacuum. Note 
that we do not take into account higher loop corrections. In the SM, the constraints 
become significantly weaker if higher loop corrections are taken into account [33,34]. 

3.2.2 Experimental constraints 

Constraints from the direct search for W and Z' 

Since the production cross sections of the extra vector bosons are relatively large, this 
model is constrained from current results of the exotic resonance searches of various 
decay channels at the y/s = 8 TeV LHC. We take account of the following constraints: 
W —)■ Ip searches [35,36], Z' —>■ ££ searches [37,38], Vh resonance searches {V' —)• Vh) 
[39-42], and the diboson searches {V —)■ VV) using dijets [ 1 , 2 ], ^^jj [43], lujj [44], and 
channel [45]. The searches for the other channels do not constrain this model. We 
fill the excluded regions with blue in Figs. 5 and 6 . Among the constraints, V ^ Vh 
and V' ^ WZ —>■ iujj give severe bound, and exclude a part of the parameter regions 
in which we can explain the diboson excess reported by the ATLAS experiment. 


Constraints from the electroweak precision measurements 

The electroweak precision parameters, S, T, W and Y, defined in Ref. [46], are severely 
restricted from the electroweak precision observables. Since the interactions of W and 
Z' to the light fermions affect the low energy observables, the light W and Z' are 
severely constrained. They are calculated at the tree level in Ref. [23], 


A 


qIvIvI _^ 

gf{vj + vi)^ + g^vf m^r, 


f =0, 



(50) 

(51) 
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w 

V 


= 4^2 _ ^1 _ 

^gfvf + vigf(vf + vir + g^,vf 

= 4^2 _ ^1 _ 

^gfvf + vigf(vf + vir + gX 



(52) 

(53) 


We use them to find the constraint on the parameter space. The parameter regions 
constrained at 95% C.L. are filled with red in Figs. 5 and 6 . The small v^, the small r 
as well as the light W regions are constrained. 

We also consider constraints from flavor physics, and find they are weaker than the 
constraints from S and/or the current LHC bound. For example, K^-K^ mixing in this 
model is almost the same as that in the SM. This is because the contributions from W is 
sufficiently suppressed due to the suppression of the couplings to the SM fermions, and 
also the the modification of the W couplings to the fermions is very small, 

Therefore, we do not show the constraints in the figures. 


3.2.3 Summary of the constraints 

In Fig. 5, we plot all the constraints in r-v^ planes with three different parameter sets, 
(m^/, ruA, i^f) = (2 TeV, 2 TeV, 1.00), (2 TeV, 2 TeV, 0.99), and (2 TeV, 1 TeV, 
1.00). Here we take all the heavy scalar bosons are degenerate, = mn' = tuh- 

The colored regions are excluded or constrained, and the white regions are allowed from 
all constraints. The gray regions surrounded by the black dotted lines represent the 
perturbativity condition of the Higgs quartic couplings for A = 10 (darker) and 100 TeV 
(lighter), and the yellow regions are that of gi. The bounded below condition excludes 
the cyan region. The global minimum vacuum condition excludes the green region. The 
magenta regions are excluded by the stability condition for A = 10 (darker) and 100 
TeV (lighter). The LHC results exclude blue regions, where the solid blue lines represent 
W —)■ If searches [35,36], the dashed lines represent Z' —>■ ££ searches [37,38], the dotted 
lines represent Vh resonance searches [39-42], and the dot-dashed lines represent diboson 
searches iV' —)■ VV) [1,2,43-45]. The regions filled with the red color are excluded by 
the electroweak precision measurements. No physical solutions are found in the black 
region, namely the gauge couplings and/or the VEVs become complex numbers there. 

By comparing all the panels, we find the experimental bounds (the LHC and the 
electroweak precision measurements) are almost insensitive to the heavy Higgs bosons. 
The theoretical bounds are very sensitive to the kf, as we can see from the panels (a) 
and (b). This is because the deviation of Kp from 1 leads to the large A 3 (cf. A 3 ~ 0.13 
(0.78) at Kp = 1.00 (0.99)) (see Eq. (35)), which changes the regions excluded by the 
bounded below condition (see Eq. (48)) and stability condition (see Eq. (49)). We find 
that perturbativity condition is weaker in the panel (c) (m^ = 1 TeV) compared with 
the panel (a) {rriA = 2 TeV). This is because the lighter leads to the smaller quartic 
couplings. 

Figure 6 shows the constraints in m^/-r plane with the same color notation in Fig. 5. 
We take the universal masses mz’ = tua = itih' = '^h, and three different choices for ^3 
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mz-= = 2 TeV, Kp= 1.00 



mz=2 TeV, rnp^^„yHigg, = 1 TeV, Kp= 1.00 



(c) 

Figure 5. The theoretical and experimental constraints in r-vs planes. We take three different 
parameter choices for {mz', mA, up)- The colored regions are constrained, Gray: the per- 
turbativity conditions for the Higgs quartic couplings, Yellow: the perturbativity conditions 
for 5 i, Cyan: the bounded below condition. Green: the global minimum vacuum condition, 
Magenta: the stability condition. Bine: the LHC bounds. Red: constraints from the elec- 
troweak precision measurements. Black: no physical solutions. See also the explanations in 
the text. 
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Figure 6. The theoretical and experimental constraints in the mz'-r planes. We take three 
different parameter choices for (us, up), and the universal masses mz' = niA = 'mpi = ttih 
are taken. A color notation is the same as in Fig. 5. 
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mz' = H^ggs = 2 TcV, Kp-= 1.00 



Figure 7. Contours of the cross section a{pp —)• VF')Br(Ty' WZ) + a{pp Z')Br(Z' — 
WW) at the LHC ^/s = 8TeV in fb unit. We take mz' = ttia = = mn = 2 TeV, 

Kp = 1.00. The blue regions are excluded by the current experimental bounds, and the 
regions of gi{p = mz') > dvr are filled with yellow. 


and Kp, {vs, Kp) =(200 GeV, 1.00), (180 GeV, 1.00), and (180 GeV, 0.99). The pertnrba- 
tivity condition for the Higgs qnartic conplings gives severe bonnds in the heavy Higgs 
mass region. This is becanse the large mass differences in the GP-even Higgs mass 
spectra reqnire the large Higgs qnartic conplings, and they become non-pertnrbative 
eventnally. We also find in the panel (c) that when Kp deviates from 1, the bounded be¬ 
low condition gives stringent constraint. This is because that the Higgs qnartic couplings 
\i are sensitive to the small deviation of Kp, see Appendix A. 

3.3 Current status: 8 TeV analyses 

In this subsection, we focus on the cross section times branching ratios of V at the LHG 
8 TeV. Inspired by the recent ATLAS diboson excess [1], we concentrate on the case 
that the masses of the extra vector bosons are around 2 TeV. 

Since the separation of the WW, WZ, and ZZ channels are not good and there 
are significant overlap among them [1], we investigate the total cross section a{pp —)■ 
W')Bi{W' WZ) + a{pp Z')Bi{Z' WW) at = 8 TeV in Fig. 7. In this 
figure, mz' = tka = mp' = rnp = 2 TeV, Kp = 1.00 are taken. The blue regions are 
excluded by the current experimental bounds discussed in Sec. 3.2, and the regions of 
gi{p = mz') > dvr is hlled with yellow. Here we do not show the constraints from the 
Higgs sector. 

We have to estimate the cross section value, a{pp —)■ V' —)■ VV), required for the 
explanation of the diboson excess. We hnd a = 6 fb with large error when we use the 
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cr(pp ^ W'^ WZ)+o-(pp ^ Z'^ WW)[fb] 


o-(pp -> W’^ WZ)+o-(pp ^ Z'-> WW)[fb] 


o-(pp ^ W’^ WZ)+o-(pp ^ Z’^ WW)[fb] 



Figure 8. The cross sections of W' and Z' with two gauge bosons in their final states. We take 
niz' = 1.9 TeV (the left panel), 2.0 TeV (the middle panel), and 2.1 TeV (the right panel). 
The color filled regions are excluded or constrained. The blue regions are excluded by the 
experimental bounds. The regions of gi(/U = mz') > dvr are represented by yellow. The cyan 
regions are excluded by the bounded below condition. The green regions are constrained from 
the perturbativity and stability conditions where their cut off scale is 100 (10) TeV in the lighter 
green (darker green) regions. Below the thick orange line the cross section a{pp —?■ W —)• VV) 
is consistent with the ATLAS result. After taking account of the iL-factor, the boundaries of 
the regions change into the dashed orange lines, and the regions below the dashed blue lines 
are excluded by the LHC bounds. 


event numbers between 1.85 and 2.15 TeV bins, the estimated background, and the 
efficiency given in Ref. [1]. Hereafter we require a = 6 fb for the explanation of the 
diboson excess. 

This cross section value can be achieved in the regions where r is much smaller than 
1. This is because the production cross section is enhanced by r“^. Hence, we focus on 
r 1 regions in the rest of this paper. 

Note that the strongest LHC constraint to the parameter regions where the ATLAS 
diboson excess can be explained comes from the hadronic channel of a{pp —?• V' —)■ 
Vh) [42], e.g. a{pp —)■ V' —)■ Vh) < 7 fb for mz' = 2 TeV. This implies a{pp — )■ V' — )■ 
VV) < 7 fb as we discussed in Sec. 3.1. If the coupling ratio kz deviates from 1, this 
LHC bound becomes weaker, because the relation between a{pp —)• V' — )• VV) and 
a{pp —)■ V' — )■ Vh) are modihed. However, in this model kz is severely constrained close 
to 1 in the rriA S> mh regime (see Appendix A). 

We show a{pp —)■ V' —)■ VV) in Fig. 8, for mz' = 1.9, 2.0, and 2.1 TeV. We take 
Kp = 1.00, and all heavy scalar masses to be the same as m^/, mz' = mA = ’mn' = rnu- 
The color hlled regions are excluded or constrained. The color notation is given in the 
caption. We hnd that the cross section is sensitive to the mass. For example, by changing 
mz' from 1.9 to 2 TeV (5 % mass difference), the cross section is decreased by about 
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0-(pp ^ W’^ Wh)[fb] 


cr(pp -» W’^ Wh)[fb] 


o-(pp ^ W’^ Wh)[fb] 


220 


> 200 
O 


Figure 9. The cross sections of the diboson channels with h. We take mw' — 1-9 TeV (the 
left panel), 2.0 TeV (the middle panel), and 2.1 TeV (the right panel). The colored regions 
are experimental and theoretical constraints. 


40%. This is because the PDF rapidly changes in the heavier mass regions (see Fig. 2). 

Below the thick orange line the cross section a{pp VV) is consistent with 

the ATLAS result. Here we apply the event selection efficiencies for the extended gauge 
model (cf. 10-16 % at mjj = 2 TeV ) [1]. Note that a part of these regions are 
constrained from the stability condition at A = 100 TeV filled with light green. However, 
once we take account of the higher order correction, these constraint would be weaker as 
we discussed in Sec. 3.2.1. In the parameter regions shown in the figure, a{pp —)■ W —)■ 
WZ)/a{pp —)■ Z' —H^H^) ~ 2. This is because the custodial symmetry is enhanced in 
the small r regime. 

In the figure we show the leading order (LO) production cross sections. Next-to- 
leading order and next-to-leading logarithmic (NLO-I-NLL) corrections to the production 
cross sections of W and Z' are evaluated in Ref [47,48], and the iF-factor (ct/ctlo) is 
about 1.3. This means that once we consider the QCD corrections, the production cross 
sections in the hgures should be scaled by about 30 %, and the LHC bounds become 
severer, while the theoretical constraints do not change. After taking account of the 
A-factor, the LHC diboson excess can be explained for the regions below dashed orange 
lines, and the regions below the dashed blue lines are excluded by the LHC bounds. 

Figure 9 shows a{pp —?• W —?• Wh).^‘^ The parameter choices and the color notations 
are the same as Fig. 8. We find that a{pp —)■ hF' —)■ Wh)/a{pp —)■ Z' — Zh) ~ 2 due to 
the enhancement of the custodial symmetry. We find a{pp —)■ W Wh) x20 fb“^ x 
Br(h —)■ hh) Br(fF —>■ eu + pu) ~ 9 events with a ~ 4 fb in the regions where the ATLAS 
diboson excess can be explained. This is consistent with the excess of the event for the 
1.8-1.9 TeV bins at the CMS with a local signihcance of 2.2a for W —)■ Wh —)■ li>hh 
search [41], although the detail of the event selection has not been reported. 

^"^Detailed studies of this process are found in Refs. [49,50]. 
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Figure 10. rtot(lF'), {gw’wzlg^z) ^ {^w'/rnwY, and gw'ff/g^ff for niz' = 2 TeV. 
They are insensitive to small difference of mz'- The color notations are the same as in Fig. 8. 


We show rtot(W^')) i^z, {gw'wz/g^wz) X {mw'/mwY and gwff/g^ff in Fig. 10. 
We fix mz' = 2 TeV because they are not sensitive to the Z' mass. The choices of other 
parameters and the color notations are the same as in Fig. 8. The width of W is shown in 
the top-left panel in the figure. We find that it is narrow and is less than 1 % of its mass, 
because the W couplings to the SM particles are suppressed by powers of mw/mw' and 
the decay into the heavy scalars are suppressed kinematically. The observable related 
to the Higgs couplings kz defined in Eq. (28) is shown in the top-right panels. The 
deviation from the SM prediction is small and the model is consistent with the current 
LHC data [51,52]. Since the International Linear Collider (ILC) can measure the kz at 
1% level [53], some parameter points are within the reach of the proposed ILC. We show 
gwwz/gy^z gwff/g^ff at the lower panels in the figure. In the benchmark model 
used in Ref. [1], gwwz/g^wz = and gwff/g^ff = 1- In our model, we find, 

due to the extra suppression by and small r, gw'wz/gwwz 1^ numerically 

the same order as the benchmark model although its niw' dependence is mw/mwi (see 
Eq. (27)). On the other hand, gw'ff/g^ff is about 10% of the benchmark model. 
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4 MC simulation of W' WZ at y/s = 13 TeV 

In this section, we perform a collider simulation of pp —>■ W WZ aX yfs = 13 
TeV. To study a discovery potential, we generate both QCD dijet background and 
pp —)■ W WZ signal events. 

We generate 1.73 x 10® QCD dijet events as the dominant background by using 
Pythia 8.205 [54] with the generation cut so that the parton-parton center of mass 
energy must exceed 1 TeV and pr > 400 GeV at a/s = 13 TeV. The tree level production 
cross section is 350 pb. Our sample therefore corresponds to roughly j dtC =5 fb“^. We 
use the Tune 4C for fragmentation and hadronization [55]. We also generate 10^ signal 
events {pp —)■ W'{W'^ + W'~) —)■ WZ ) for the mass between 1800 and 3200 GeV, and 
take rtot(hT') to be 25 GeV in the simulation. Note that the total width of W is less 
than about 30 GeV in the allowed region of this model (see Fig. 10). We also generate 
the signal and background events at yfs = 8 TeV and compare them with the ATLAS 
plots [1]. 

The simple detector simulator DelphesS [56] is modihed using FastJetS [57,58] so 
that the mass drop and the grooming cuts used in the ATLAS study can be applied to the 
jets. We apply the cluster track matching algorithm of DelphesS so that information 
of tracks inside jets can be used, otherwise the default ATLAS card is used. 

Reconstruction of boosted objects using jet substructure was originally proposed in 
Refs. [59,60]. See recent developments in Refs. [61,62]. In our simulation, we closely 
follow the ATLAS analysis. The Gambridge-Aachen algorithm with R = 1.2 is used 
[63,64] for the jet clustering. Then, pn > 600 GeV, pTa > 540 GeV, (pti — Pt 2 )/{PTi + 
PT 2 ) < 0.15, |pi — P 2 I < 1-2, |pi| <2 and |? 72 | < 2 are required for the jets. In addition, 
we require Exmiss < 350 GeV, and veto events with isolated electrons and muons with 
Pt > 20 GeV. For each jet, the pair of the subjets which satishes the subjet momentum 
balance criteria ^/p > = 0.45 are selected, where 

j — . / \ A \ 

Vv = mm Pt,, , Pt,^ -■ 54 

mo 

Here, pt,, and pt,^ are the transverse momenta of subjets ji and ^ 2 , ^R{ji,j 2 ) is the 
distance between the subjets ji and j 2 , and mo is the mass of the parent jet. Then 
the constituents of the selected pair of subjets are hltered. Namely the constituents are 
clustered with the radius parameter R = 0.3, and up to the highest 3 jets are taken to 
calculate the groomed jet mass and momentum. We require |my — m/ < 13 GeV, where 
my is ruz or mw, and rrij is an invariant mass of the groomed jet. Finally, the number 
of charged-particle tracks which are associated with the jet is required to be rich < 30. 
In Fig. 11, we show the distribution in rich and mj plane for the events with mjj > 1500 
GeV where all the cuts except for rich and rrij are applied. The ATLAS signal regions 
are marked by squares. The hgure shows very good separation between the signal and 
the background events. 

To check our simulation, we compare the distributions of our a/s = 8 TeV samples 
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Figure 11. The distributions of the signal {mw' = 1800 GeV and rtot(hF0 = 25 GeV) and the 
dijet background in rich and mj plane. All the cut except for rich and rrij are applied, and we 
required rujj > 1500 GeV. The ATLAS signal regions are marked by squares. 


to the ATLAS ones. The Fig. 12 shows the reconstruction efficiency of the signal for 
various input W mass. Our result at \/s = 13 TeV is also shown. We hnd that the 
signal selection efficiency agrees with the ATLAS one. Since the jets get narrower with 
increasing the signal efficiency becomes lower for higher W' mass. 

On the other hand, we hnd that the number of the background events after all the 
selection cuts are twice as large as that of the ATLAS hnal result. The discrepancy in 
the total background events might arise from several sources. The number of charged 
tracks of QCD jets are controlled by soft physics and varies signihcantly depending on 
Monte Carlo (MC) generators and tunes of the parton shower parameters especially for 
gluon jets. The distributions are shown in the left panel of Fig. 13, under the cut of 
Fig. 1 of Ref. [1], 1.62 TeV< mjj < 1.98 TeV and 60 GeV < mj < 110 GeV together 
with the selection cuts listed above except that for rich, where rrijj is the dijet invariant 
mass. The signal distribution which is represented in the blue line agrees quite well 
with the ATLAS ones. However, the average number of charged tracks of dijet event 
is signihcantly higher than the ATLAS ones. To see the MC dependence, we also show 
the distribution of the MC sample generated by Herwig++ with default tunes by the 
black line [65], which predicts slightly small (uck) compared with PythiaS Tune C4, 
but higher than ATLAS CTIO Tune.^® 

In the right panel of Fig. 13, we also show the mj distributions for the signal and 
background. The signal distribution agrees quite well with the ATLAS MC results again. 
For the background distribution, we hnd that the number of events above mj > 50 GeV 
is smaller compared with Fig. 1 of Ref. [1]. This is because we do not generate underlying 
events together with the dijet events. 

In the experimental side, ATLAS counts the well reconstructed track inside the 
jet. The efficiencies are not implemented in our simulation. Naively speaking, the 

In the study of quark-gluon separation [66], HerwigH — h reproduces high gluon pr jet nature well. 
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Figure 12. The event selection efficiencies for W' —)• ITZ —>■ JJ simulated events generated by 
Pythia 8 as a function of the W' mass at ^/s = 8 and 13 TeV. The red band is the result of 
the ATLAS [1], where the thickness corresponds to ± 1 cr statistical and systematical errors. 
Our result of MC simulation at y/s = 8 (13) TeV is represented in the blue (green) band where 
the thickness corresponds to ± 1 u statistical error. 


efficiency is expected to be lower for the jets with high charged track multiplicity. For 
jet clustering, ATLAS selects calorimeter towers using Topocluster algorithm and does 
not use cluster-track matching which is only crudely implemented in our simulation. 
In any case, data driven approaches are adopted in Ref. [1] to estimate the selection 
efficiency, and reproducing the result precisely is beyond the scope of this paper. It is 
probably worth doing more dedicated theoretical and experimental studies on jet nature 
relevant to the boosted W and Z bosons reconstruction in future. 

Keeping crudeness of our simulation in mind, we estimate the signal efficiency at 
y/s = 13 TeV using our signal MC and detector simulation without rescaling, while the 
number of the background events obtained from our MC is rescaled by factor of 1/2 
which is needed to reproduce the ATLAS results at y/s = 8 TeV. The scaling approach 
comes from an assumption that the change in the center of mass energy from 8 TeV to 
13 TeV does not alter the structure inside jets of the same pt jets. Note that gluon jets 
are involved in the QCD background at y/s = 13 TeV but this effect is not taken into 
account. Under this assumption, the background is htted to estimate the distribution, 
and the result is shown in the left panel of Fig. 14. We found 169 events/10 fb“^ for 
1550 GeV < rrijj < 3550 GeV after rescaling. We also show the dijet invariant mass 
distributions of the signal for various input W mass in the right panel of Fig. 14. 

Figure 15 shows expected the exclusion limit at 95 % C.L. for a{pp —>■ U' —>■ VV) 
at \/s = 13 TeV."^® We calculated of the signal plus background distributions to 
the background distribution. For the signal, we generate 10^ signal events to obtain the 
result for each mw', and take the number of events in the bins imax — 3 < i < imax + 3, 

^®Prospects for the electroweak gauge boson scattering which can also probe W are discussed in, for 
example, Ref. [67]. 
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Number of charged track 


jet mass 




Figure 13. The distributions of the number of charged tracks rich (left) and the jet mass rrij 
(right) for the signal (Blue) and background events (Black and Red). We take mw' = 1800 
GeV and rtot(bF0 = 25 GeV. We simulate the background events using PythiaS (Red) and 
HERWIG++ (Black). 


where Vax is the highest signal event bin and the bin size is 50 GeV. Then, we rescale 
them by the cross section and Inminosity. The selection efficiency of the signal event is 
shown in Fig. 12. For the backgronnd, we nse onr £t and rescaled it by the Inminosity. 
Here y-axis means a{pp —)■ W —)■ IFZ) + o{j)p —)■ Z' —)■ IFhF) with the mass at myi. 
We find that the region where a{pp —)■ V' —)■ VV) is larger than 20 fb may be exclnded 
at J dtC = 10 fb~^ and ^/s = 13 TeV. 

5 Future prospects: 13 TeV analyses 

In this section, we analyze the fntnre prospects of W and Z' searches by applying result 
in Sec. 4. 

We investigate the prospects for the y/s = 13 TeV collision in the case of mz' = 2 
TeV. In Fig. 16, we show contours of the several cross sections by orange dashed lines 
in r-us planes: a{pp —)■ W')Bi{W' WZ) + a{pp —)■ Z')Bi{Z' —)■ WW), a{pp —)■ 
Wh), a{pp W')Bt{W' iu), a{pp W')Bt{W' jj) + alpp 
Z')Br(Z' —)■ jj), and a{pp — )■ Z')Bt{Z' —)■ ££). The color Filed regions are excluded and 
constrained as we discussed in Sec. 3.3. The color notations are the same as in Fig. 8. 

The expected exclusion limit at \/~s = 13 TeV provided in Fig. 15 is shown by the red 
(dashed) line for j dtC = 10 (100) fb“^. The masses of W and Z' are highly degenerate 
in the small r regime, and we assume that the signal efficiency for Z' —)■ WW event is 
equivalent to the one for W —)■ IFZ which is simulated in the previous section. Thus we 
can apply the prospect shown in Fig. 15 to a{pp —)• W WZ) + a{pp —)■ Z' —)■ WW). 
The region where the ATLAS diboson excess can be explained is within the reach of the 
LHC at f dtC = 10 fb~^ and ^/s = 13 TeV. In addition, the cross sections of the other 
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reconstructed trijj distribution at 13TeV 




Figure 14. The rrijj distributions for the background (left) and the signal (right) at y/s = 13 
TeV. The background distribution is for 10 fb“^, and the signal distributions are normalized 
to be 1 for various input W mass. 


Expected limit at 13 TeV 



my' [GeV] 

Figure 15. Expected limit on a{pp —)• V' —)• VV) at ^/s = 13 TeV. We take Ftot(V') = 25 GeV. 
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Figure 16. The prospects for the ^/s = 13 TeV collision. We take mz' = tua = mn' = fnn = 
2000 GeV, Kp = 1.00. The color notation is the same as in Fig. 8. Below the thick orange 
line, the cross section a{pp —)■ W —)• VV) is enough to explain the diboson excess at ATLAS. 
Below the (dashed) red line is expected to be excluded with f dtC = 10 (100) fb“^. 


28 










































channels are also large, so that the spin-1 resonances could be probed in the channels 
as well. 

We plot the a/s = 13 TeV cross sections as a function of spin-1 resonance mass 
in Fig. 17. We take the universal mass m^/ = uia = niH' = '^H) = 200 GeV, 

and Kp = 1.00 in Fig. 17. The color notation is the same as in Fig. 8. The mass of 
W is highly degenerate with the mass of Z' in this hgure. We hnd that the spin-1 
resonances lighter than 2.1 (2.6) TeV can be excluded by the diboson search at the 
LHC with f dtC = 10 (100) fb“^. Figure 18 shows the same cross sections for another 
parameter set = 150 GeV and Kp = 0.99). With the consideration of the bounded 
below condition, we hnd that the spin-1 resonances up to masses of 2.3 (2.6) TeV can 
be excluded with J dtC = 10 (100) fb“h 

Before closing this section, we briehy comment on the prospect for the large r region. 
Under the several theoretical bounds, this model has two allowed regions, namely small 
r and large r regions (see Fig. 5). The cross section of pp ^ Z' ^ ll is enhanced by 
factor in large r regime (see Eq. (72)). ATLAS study found around 3 TeV with the 
cross section around 0.01 fb is accessible with the high luminosity LHC at ^/s = 14 TeV 
with integrated luminosity 3000 fb“^ [68]. Thus we can expect that the cross section for 
the large r region in the model is accessible at the LHC Run-2 as we can see in Fig. 16. 
However, the r regime would be constrained from the electroweak precision parameters 
at one-loop level due to the custodial symmetry breaking. 

6 Conclusion 

Motivated by the ATLAS diboson excess around 2 TeV, we have investigated the phe¬ 
nomenology of the spin-1 resonances iW' and Z') in the partially composite standard 
model. In this model, W and Z' couple to the SM fermions weakly through the mixing 
to the elementary gauge bosons. We hnd that the main decay modes of the resonances 
are V —)■ VV and V —)■ Vh, and the width is narrow enough so that the ATLAS 
diboson excess can be explained. The couplings of the spin-1 resonances with the SM 
sector can be controlled by the ratio of the Higgs VEVs so that the ATLAS diboson 
excess can be explained. 

We have explored not only the current bounds from the LHC and the precision mea¬ 
surements but also the theoretical constraints, i.e. perturbativity condition, bounded 
below condition, global minimum vacuum condition, and stability condition of the scalar 
potential. The parameter regions where the diboson excess at the ATLAS can be ex¬ 
plained are still allowed after including those constraints. 

In order to investigate future prospects of the spin-1 resonance search, we have 
performed the simulation at a/s = 13 TeV LHC, and estimated model independent 
exclusion limit for a{pp —)■ W —)■ VV JJ) shown in Fig. 15. Applying our simulation 
result, we hnd that the parameter regions consistent with the ATLAS diboson excess 
will be excluded at J dtC = 10 fb“^ and a/s = 13 TeV. 
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o-(pp ^ W'^ WZ)+0-(pp ^ Z'^ WW)[fb] 



mz[GeV] 


cr(pp ^ W’^ Wh)[fb] 



mz[GeV] 


(r(pp ^ Z’^ l()m 



mz[GeV] 


Figure 17. The prospects of the discovery as a function of mz' and r for the LHC Run-2. We 
take U 3 = 200 GeV and Kp = 1.00. The color notation is the same as in Fig. 8 . The (dashed) 
red is the future expected exclusion limit for f dtC = 10 (100) fb“^. 



o-(pp -> Z'-> /’/’)[fb] 


cr(pp ^ W’^ Wh)[fb] 


mz[GeV] 


mz[GeV] 


Figure 18. The prospects of the discovery as a function of mz' and r for the LHC Run-2. We 
take V 3 = 150 GeV and Kp = 0.99. Below the (dashed) red line is expected to be excluded 
with f dtC = 10 (100) fb“^. 


30 





























Finally, we have investigated future prospects of diboson resonance search in our 
model. The spin-1 resonances up to a mass of 2.6 TeV can be probed at y/s = 13 TeV 
and J dtC = 100 fb“^. 
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A Viable range of the coupling ratio kz 


In this appendix, we show kz is very restricted to be one. We hrst replace the Higgs 
quartic couplings A* by the other parameters Eq. (34). Here we take uia = ttih' = rnn 
for simplicity, and work in ruA S> mt regime. For r 1 regime, we hnd 


and for r 3> 1 regime. 


Ai ~0, 

A 2 — K,z)t 
, ml n m\ . 2 N 

Ai2 ^0, 

•^23 —/2(«F, «^z), 

A 31 ^ 0 , 


Al 

A2 ~ 0 , 

, ml o m\ 2 N 

^3=-^'‘f+^(1-Ff), 

Ai2 ^0, 

A23 -0, 

A3I —/2(kf, Kz)i 


where 


ml 


/i(kf, k,z) —^ ~ (1 “ i^f) 


2 {m 


ml 


+ (l-«:z)- 


,2 ’ 


(55) 

(56) 

(57) 

(58) 

(59) 

(60) 


(61) 

(62) 

(63) 

(64) 

(65) 

( 66 ) 
(67) 


( 68 ) 
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mz=2000GeV, V3=200GeV 
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r 

Figure 19. The viable range of the coupling ratio kz, where we take us = 200 GeV, mz' = 
vriA = rnu' = ran = 2 TeV, and Kp = 1-00. The regions where Aj(^ = 10 (100) TeV) > (dvr)^ 
are filled with green (light green). 



ml 


/2(k_f, ^z) —~ + (1 

/jiZ 


Kp) 


— 3^3 m\ — m\ 


+ (1 — ^z) 


m\-ml 


(69) 


In Eqs. (68) and (69), we make an expansion around ~ 1 and kz — 1. Since both 
coefficients of (1 — Kp) and (1 — kz) are large enough in 3> v regime, even the small 
deviations oi Kp and kz from 1 change the Higgs quartic couplings Aj drastically. 

Let us consider the case of Kp = 1. The largest Higgs quartic coupling is A 23 (A 31 ) 
in the r -C 1 (r i:g> 1) regime. We further demand A 23 (A 31 ) < (dvr)^, namely 


Kz > 


(167rV -m^) A _ ^ 
{m\ - ml) V 


(70) 


Due to the running effects, these quartic couplings can be even larger at the high scale, 
so that the lower bound on kz is actually severer than above estimation. 

In Fig. 19, we show that the viable range of the coupling ratio kz, where we take 
^3 = 200 GeV, mz' = m^ = mn' = mp = 2 TeV, and Kp = 1.00. The regions where 
Aj(/i = 10 (100) TeV) > (dvr)^ are hlled with green (light green), and there is no physical 
solution in the gray region. The maximal value of kz is achieved at a boundary of the 
gray region. The white area represents the allowed region of parameter space. Thus we 
hnd Kz is severely constrained close to its maximized value for ^ mh- 
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B Production cross sections of W' and Z' 

The leading order production cross sections of W and Z' are given as follows. 


a{pp^W'^X) ~ 


a{pp —)■ Z'X) ~ 


71 1 

12s s^ r2 


i-d 


dx 




X 


X fu{x,mw')fd 


m 


W' 


sx 


,mw' + fd{x,mw')fu 


m 


W' 


sx 


+fd{x,mw')fu 


m 


W' 


sx 


,mw' + fu{x,mw')fd 


m 


W' 


sx 


, miy/ 


, mw' 


T {u c, d ’ 


(71) 


71 1 

6s s^ r2 


1-^ 


X 




s2 


J- Tf+^ZQf +(^"ZQf 


s2 


r dx 

' I s X 


m 


X ( fu{x,mz')fu ( '—,mz' ) + fu{x,mz')f, 


sx 


mz' 

sx 


,mz' 


+ {u ^ d, s, c, b)], 


(72) 


where s is a square of the center of mass energy of the pp collider and fq{x, Q) is the 
parton distributions inside the p at the factorization scale Q for quark flavor q. Note 
that Eq. (71) is a sum of production cross section of W'^ and W'~. 


C Renormalization group equations 

We derive the one-loop [3 functions for this model [69,70], and obtain. 


/3go 

— ^9oy 

(73) 

f^gi 

= -791 

(74) 

/^92 

= 7gl 

(75) 

^gs 

= -791 

(76) 

/^Ai 

= 2A\\ + 2 A ^2 + 2 A 31 - 9Xi{gl + g^) + -g^ + -gf + -g^gj, 

(77) 

/^A 2 

= 24 A 2 -I- 2A^2 + 2 A 23 ~ 3'^2(3fi'i -I- gl) + -(?! -I- -g2 + -^9i92i 

(78) 

CO 

< 

= 24 A 3 2 A 23 + 2 A 31 — 3X3{3gl gl) -gQ + -g2 + -^9o92 



+ 2 X 3 ( 39 '^ + 3yl + yl) - ^yf - ^yt - ^yt, 

(79) 
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/^Ai2 — 4A^2 + 12Ai2(Ai + A 2 ) + 4 A 23 A 31 —-Ai2(35fg + + 5f2) + (80) 

/^A23 = 4 A 23 + 12A23(A2 + A 3 ) + 4 A 12 A 31 — -A23(35'o + + 2gl) + -g^ 

+ ^23{^yt + + l/r)) (81) 

/^Aai = 4 A 31 + 12 A 3 i(Ai + A 3 ) + 4 A 12 A 23 — 2 '^ 3 l( 6 (y'o + ^9i + fl' 2 ) + ^fl'o 

+A3i(3|/4^ + 3yl + 1 /^), (82) 

Pyt = yt - ^9l - 8gs + \yl + \yl + \yl^ , (83) 

l^y, = yh(y-\9l-^9l-^9l + \yl + \yl + \y‘^^ , (84) 

Pyr = yr (y-\9l - ^9l + ^y'i + \yl + \yi^ , (85) 


where (3 functions are defined in the following notation, 

d(cOUphng) _ /^(coupling) 

dln/i (47r)2 


( 86 ) 
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